The Density Matrix Renormalization Group and the Shell Model 
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We summarize our recent efforts to develop the Density Matrix Renormalization Group (DMRG) 
method into a practical truncation strategy for large-scale nuclear shell model calculations. Following 
an overview of the essential features of the DMRG, we discuss the changes we have implemented for 
its use in nuclei. In particular, we have found it useful to develop an angular-momentum conserving 
variant of the method (the JDMRG). We then summarize the principal results we have obtained to 
date, first reporting test results for ''^Cr and then more recent test results for ^^Ni. In both cases 
we consider nucleons limited to the 2p-lf shell. While both calculations produce a high level of 
agreement with the exact shell model results, the fraction of the complete space required to achieve 
this high level of agreement is found to go down rapidly as the size of the full space grows. 

PACS numbers: 21.60.Cs, 05.10.Cc 



I. INTRODUCTION 

In the traditional nuclear shell model, the low-energy 
structure of a nucleus is described by diagonalizing the 
effective nuclear hamiltonian in an active space consisting 
of at most a few major shells outside/inside an assumed 
doubly-magic core. Despite the enormous truncation in- 
herent in this approach, it can still only be applied in very 
limited nuclear regimes. For heavy nuclei or nuclei far 
from closed shells, further truncation of the shell-model 
space to a manageable size is required. 

In this work, we discuss the use of the Density Matrix 
Renormalization Group (DMRG) as a truncation strat- 
egy for the nuclear shell model. Originally developed for 
the treatment of low-dimensional quantum lattices 
the DMRG method was subsequently extended with im- 
pressive success 0, 0, to finite Fermi systems, sug- 
gesting its possible usefulness for the nuclear shell model. 

The DMRG method involves a systematic inclusion of 
the degrees of freedom of the problem. When treating 
quantum lattices, real-space sites are added iteratively. 
In finite Fermi systems, these sites are replaced by single- 
particle levels. At each stage, the group of sites that 
has been treated (referred to as a block) is enlarged to 
include an additional site or level. This enlarged block 
is then coupled to the rest of the system (the medium) 
giving rise to the superblock. For a given eigenstate of 
the superblock (often the ground state) or perhaps for a 
group of important eigenstates, the reduced density ma- 
trix of the enlarged block in the presence of the medium is 
constructed and diagonalized and those states with the 
largest eigenvalues are retained. This method of trun- 
cation is guaranteed to be optimal in the sense that it 
maximizes the overlap of the truncated wave function 
with the superblock wave function prior to truncation. 

This process of systematically growing the system and 
determining the optimal structure within that enlarged 
block is carried out iteratively, by sweeping back and 
forth through the sites, at each stage using the results 
from the previous sweep to define the medium. In this 



way, the process iteratively updates the information on 
each block until convergence from one sweep to the next 
is achieved. Finally, the calculations are carried out as a 
function of the number of states retained in each block, 
until the changes are acceptably small. 

The traditional DMRG method works in a product 
space, whereby the enlarged block is obtained as a prod- 
uct of states in the block and the added site and the su- 
perblock is obtained as a product of states in the enlarged 
block and the medium. In the nuclear context, this means 
working in the m-scheme. This method was applied in 
nuclear physics by Papenbrock and Dean [5| and shown 
to provide an accurate description of the properties of 
'^^Si, but a much poorer description of ^^Ni where the 
converged solution was still energetically quite far from 
the exact ground state. 

A limitation of the traditional algorithm is that it 
does not preserve symmetries throughout the iterative 
process. Since the density matrix procedure involves a 
truncation at each iterative stage, there is the potential 
to lose these symmetries and the associated correlations. 
On this basis, we proposed Q the adoption of a strategy 
in which angular momentum is preserved throughout the 
DMRG process. This method, called the JDMRG, was 
applied in nuclear physics for the first time in the con- 
text of the Gamow Shell Model 0] and then subsequently 
developed as a means of approximating the traditional 
nuclear shell model, with a preliminary application first 
reported [3] for '^^Cr. More recently, the algorithm was 
computationally improved and applied both to '^^Cr and 
to the ground state of ^^Ni Q. In this work, we extend 
the results of ref. Q to include excited states of ^^Ni. 

An outline of the presentation is as follows. We be- 
gin in Section II with a brief overview of the traditional 
DMRG method including a discussion of the changes re- 
quired for exact angular-momentum conservation. In 
Section III, we report our results for ""^Cr, which were 
first presented in j^] , and then report our most recent re- 
sults for ^^Ni, including those for excited states. Finally, 
in Section IV we summarize the principal conclusions of 
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FIG. 1; Schematic illustration of the DMRG growth proce- 
dure. A block B consisting of sites 1 and 2 is enlarged to 
include site 3, forming B' . The medium M consists of all of 
the remaining sites, 4 — 8. 



this work and outline some directions for future investi- 
gation. 



II. AN OVERVIEW OF THE DMRG METHOD 
A. The truncation strategy 

The DMRG method is based on an iterative inclusion 
of the degrees of freedom of the problem, represented 
as sites on a lattice. This is illustrated schematically in 
Figure 1 for a system with 8 ordered sites. 

Assume that we have treated a group of sites, referred 
to as the block and denoted B, and that we have retained 
a total number of (optimal) states m within that block. 
We now wish to add to this block the next site (r) with 
I states, thereby producing an enlarged block B' . For 
the moment, we will assume a product (or m-scheme) 
description, so that the enlarged block has m x I states. 



\hj >B'= \i >B \j >r 



i = l,m , j ^ 1,1 . (1) 



We wish to retain the optimal m states for the enlarged 
block. How do we choose them? 

In the DMRG method, we consider the enlarged block 
in the presence of a medium M that reflects all of the 
other sites of the system, referring to the full system as 
the superblock (SB). Assuming that the medium is like- 
wise described by its optimum m states, the m x I x m 
states of the superblock can be expressed as 



k >SB= \i,j >B' \k > 



M 



(2) 



We then diagonalize the full hamiltonian of the system 
in the superblock, isolating on its ground state. 



\GS >SB= ^ 'i'ijk\i,j,k >SB 



(3) 



If we then construct the reduced density matrix of the 
enlarged block in the ground state, 



Pij,i'i' 



k 



(4) 



diagonalize it and retain the m eigenstates with the 
largest eigenvalues we are guaranteed to have the m most 
important (or optimal) states of the enlarged block in the 
ground state ([3]) of the superblock. 

It is straightforward to target a group of states of the 
system, and not just the ground state, by constructing 
a mixed density matrix containing information on all of 
them. 

Once the optimal m states are chosen, we renormalize 
all required operators of the problem to the truncated 
space and store this information. This includes all sub- 
operators of the hamiltonian. 



a|, ajaj 



a] a], 



aja]afc, a|a]afea;. 



h.c. 



Having this information for the block and the additional 
level or site enables us to calculate all such matrix ele- 
ments for the enlarged block as needed in the iterative 
growth procedure. 



B. Stages of the DMRG method 

With the above remarks as background, the DMRG 
procedure involves the following stages. 



1. The warmup stage 

In the warmup stage we make an initial guess on the 
optimal TO states for each block. This choice is important 
in determining how rapidly the iterative procedure con- 
verges. In our treatment, we do this by growing blocks 
from each side of the chain gradually, using those orbits 
that have already been treated on the other end as the 
medium. This is illustrated schematically in Figure 2 for 
two successive steps of the warmup procedure. 



2. Ttie Sweep stage 

In this stage of the process, schematically illustrated 
in Figure 3, we gradually sweep back and forth through 
the sites of the lattice, at each step of the sweep using 
for the medium the results either from the warmup phase 
(during the first sweep) or from the previous sweep stage. 
The sweep process is done over and over until convergence 
is achieved from one sweep to the next. 



3. As a function of m 

The above calculations are done for a given choice of 
TO. The calculations are then done as a function of to. 
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FIG. 2: Schematic illustration of two successive steps in the 
warinup procedure, (a) A block involving the site 8 is enlarged 
to include the next site 7, with the medium being a block 
involving the optimum states from the sites 1 and 2 obtained 
in the previous warmup step, (b) A block involving the sites 1 
and 2 is enlarged to include the next site 3, with the medium 
being a block involving the optimum states from the sites 8 
and 7 obtained in the warmup step of part (a). 



FIG. 3: Schematic illustration of two successive steps in the 
sweep procedure, (a) A block involving the site 8 is enlarged 
to include the next site 7, with the medium being a block 
involving the optimum states from the sites 1—6 obtained 
previously, (b) A block involving the sites 8 and 7 is enlarged 
to include the next site 6, with the medium being a block 
involving the optimum states from the sites 1 — 5 obtained 
previously. 



until the changes with increasing m are acceptably small. 

C. The JDMRG approach 

As noted earlier, most DMRG approaches violate sym- 
metries. In nuclei, for example, they typically work in 
the m-scheme. Such a procedure is potentially problem- 
atic when imposing truncation, however, as it is diffi- 
cult to ensure that the states retained contain all the 
components required by the Clebsch Gordan series to 
build states of good angular momentum. For this rea- 
son, we have chosen to develop an angular-momcntum- 
conserving variant of the DMRG method in which angu- 
lar momentum is preserved throughout the growth, trun- 
cation and renormalization stages. The most significant 
difference between this (the JDMRG) approach and the 
more traditional DMRG approach is that now we must 
calculate and store the reduced matrix elements of all 
sub-operators of the hamiltonian, 

4, [aldj^, [ala]f, ([ala]]^a,)^ 
([ala]f [afca,]^)° + h.c. 
This can be done using standard Racah techniques. 



D. A three-block JDMRG strategy 

In the nuclear shell-model calculations we will report 

here, wc adopt a three-block strategy for the enlarge- 
ment and truncation process, schematically summarized 
in Figure 4. 
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FIG. 4: Schematic illustration of the three-block DMRG 
growth procedure for a system with neutron and proton levels. 

We begin by choosing our order of sites so that neutron 
and proton orbitals sit on opposite ends of the chain. Wc 
then gradually grow blocks of each type of particle only, 
namely we grow neutron blocks and proton blocks but 
no mixed blocks. Lastly, in the sweep stage wc go back 
and forth through the orbits of a given type of particle 
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only. As can be seen from the figure, the medium in 
this approach involves two components. If, for example, 
we are enlarging a proton block, the full medium (M) 
involves all of the remaining proton levels (Mi) and all 
of the neutron levels (M2). 

III. CALCULATIONS 

We have carried out test calculations of the JDMRG 
method described above on the nuclei ^*Cr and ^^Ni. 
We assume that these nuclei can be described in terms of 
valence neutrons and valence protons outside a doubly- 
magic ^°Ca core. In both, we use the order of single- 
particle levels shown in Figure 5. 

We report the results for these two applications in the 
following subsections. 

A. Results for "^Cr 

We begin with our results for ^^Cr, for which there are 
four neutrons and four protons outside a '^'^Ca core. In 
these calculations we assumed a KB3 effective interac- 
tion between valence nucleons and compared our results 
with those obtained for the same hamiltonian in ref. [13] ■ 
The size of the full shell-model space in this case involves 
1,963,461 states, of which 41,355 are 0+ states, 182,421 
are 2+ states, 246,979 are 4+ states, etc. 

Our results for the ground state are presented in Ta- 
ble I. The exact calculation produces a ground state en- 
ergy of -32.953 MeV. The DMRG calculation converges 
smoothly to this result as m is increased, but requires the 
inclusion of a substantial fraction of the full space to ob- 
tain a high level of accuracy. With of order 25% of the 
full 0+ space, we are able to achieve accuracy to only a 
few keV. To achieve an accuracy of better than 50 keV 
still requires, however, over 20 % of the full 0"*" space. 

In Figure 5, we show how the results converge as a 
function of the number of sweeps. As can be readily 
seen, after the first sweep we are extremely close to the 
final converged result. This is a direct consequence of our 
use of a warmup procedure that incorporates (in step 0) 
a significant part of the correlations. 



TABLE 1: Calculated ground-state energies in MeV as a func- 
tion of m for ^**Cr. The maximum dimension encountered in 
the sweep process is also given. 



m 


Eos 


Adax Dim 


40 


-32.698 


1,985 


60 


-32.763 


2,859 


80 


-32.788 


3,765 


100 


-32.817 


4,494 


120 


-32.840 


6,367 


140 


-32.890 


8,217 


160 


-32.902 


9,978 


180 


-32.944 


11,062 


200 


-32.947 


12,076 


Exact 


-32.953 


41,355 



In Table II, we present results for the lowest excited 
states, obtained using the blocks obtained at the ground 
state minimum. Here too convergence to the exact re- 
sults is achieved, if a sufficiently large fraction of the full 
shell-model space is retained. This is true despite the 
fact that the density matrix truncation procedure imple- 
mented targeted the ground state only. 



B. 5®iVi 

Next we turn to ^^Ni, for which the size of the full 
shell-model space in the m-scheme contains 1,087,455,228 
states, of which 15,443,685 have J"" = 0+ , 71,109,189 
have ,r = 2+ and 105,537,723 have ,r = 4+ , in all cases 
much larger than in "^^Cr. 

In our earlier published work on ^^Ni, we used the KB3 
interaction for ^^Ni as well and compared the ground 
state energy with the results from ref. [ll|. We only 
considered the ground state in those calculations, since 
those were the only exact results that had been reported 
in the literature. The reason for that is that KB3 is 
known to produce poor agreement with the experimen- 
tal spectra for nuclei in this region of the 2p — If shell. 
Much better agreement can be obtained with the im- 
proved GXPFlA interaction [l^. Thus, we have redone 
our test calculations of ^^Ni with this interaction, now 
comparing with the exact results not only for the ground 
state but for low- lying excited states as well [l3| ■ These 
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FIG. 5: Order of single-particle levels used in the DMRG 
growth algorithm for both ^^Cr and Ni. 



TABLE II: Results for the excitation energies in MeV of the 
lowest 2+ and 4+ states in **Cr from the JDMRG calcula- 
tions described in the text. The dimensions of the associated 
hamiltonian matrices are given in parentheses. 



m 






140 
160 
180 


0.873 ( 9,191) 
0.860 ( 12,038) 
0.855 ( 15,148) 


2.022 ( 12,442) 
1.996 ( 16,553) 
1.989 ( 21,628) 


Exact 


0.806 (182,421) 


1.823 (246,979) 
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are the largest test calculations we have carried out to 
date using the JDMRG method. 

Ground state of 48Cr with sweeps 
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FIG. 6: Convergence of the results for ^®Cr with the number 
of sweeps. 



TABLE IV: Results for the excitation energies in AleV of the 
lowest 2+ and 4+ states in '^'^Ni from the JDMRG calcula- 
tions described in the text. The dimensions of the associated 
hamiltonian matrices are shown in parentheses. 



m 


2+ 


4^ 


60 
80 
100 


2.970 ( 296,633) 
2.944 ( 345,213) 
2.942 ( 423,265) 


4.137 ( 445,898) 
4.123 ( 556,572) 
4.090 ( 701,502) 


Exact 


2.600 (71,109,189) 


3.688 (105,537,723) 



be possible to improve the quality of the description of 
low-lying excited states by targeting both the ground 
state and the first excited state in the density matrix 
truncation procedure, hopefully without losing too much 
accuracy in our reproduction of the ground state. 

IV. SUMMARY AND OUTLOOK 



The results for the ground state energy as a function 
of m are shown in Table III. Here we are able to achieve 
roughly 60 keV accuracy with just 0.5% of the full 0+ 
space. This is a significantly lower fraction of the full 
space than was required in *^Cr to achieve the same level 
of accuracy. 

This is perhaps the most encouraging result of our work 
to date. It suggests that the fraction of the space required 
to achieve a high level of accuracy with the JDMRG 
method goes down rapidly as the size of the space prob- 
lem increases. If this is confirmed with extension to even 
larger problems it would bode very well for the future 
usefulness of the JDMRG method as a practical trunca- 
tion approach for large-scale shell-model studies. 

As noted above, exact results also exist with the 
GXPFIA interaction for excited states. In Table IV, we 
present our results for the lowest 2+ and 4+ states, again 
in comparison with the exact results. These results were 
obtained by diagonalizing the associated hamiltonian ma- 
trices that derive with the blocks obtained at the ground 
state minimum. Here too the results are getting better 
with TO, albeit slowly. Agreement with the exact results 
for m = 100 is still not as good as we would like. It might 



TABLE III: Results for the ground-state energy of ^'^Ni in 
MeV from the JDMRG calculations described in the text. 
Max Dim refers to the maximum dimension of the superblock 
hamiltonian matrix. 



m 


Egs 


Max Dim 


60 


-205.603 


64,397 


80 


-205.632 


74,677 


100 


-205.643 


87,633 


120 


-205.652 


106,383 


Exact 


-205.709 


15,443,685 



In this talk, we have summarized the current status 
of our efforts to build the Density Matrix Renormaliza- 
tion Group Method method into a practical truncation 
strategy for large-scale shell-model calculations of atomic 
nuclei. Following an overview of the essential features of 
the method, we discussed the changes we had to imple- 
ment for its application to nuclei. Most importantly, we 
found it useful to develop an angular-momentum con- 
serving version of the method, the JDMRG. We then 
summarized the results we have obtained for the nuclei 
^^Cr and ^^Ni, in both cases comparing with the results 
of exact diagonalization. Both calculations were able to 
accurately reproduce the exact shell-model results. In 
the case of "'^Cr, however, this high level of accuracy re- 
quired us to retain a very large fraction of the full space. 
In contrast, we were able to achieve comparably accu- 
rate results for ^^Ni with a much smaller fraction of the 
space. The fact that the fraction of the space goes down 
rapidly with the size of the problem bodes well for the fu- 
ture usefulness of the method in even larger shell-model 
problems. 

There are several issues that we intend to explore in 
the near future. One concerns the need to determine 
through additional calculations how rapidly the fraction 
of the space required for convergence scales with the size 
of the problem. Currently we only have two data points, 
^^Cr and ^^Ni. More are needed to draw meaningful 
conclusions. 

We are also looking into the question of whether we 
can obtain comparable agreement for both the ground 
state and low-lying excited states by building a mixed 
density matrix that includes information not only on the 
ground state but also on the first excited 0+ state. 

We are also in the process of applying these methods 
to odd-mass nuclei around ^^Ni. Exact and Coupled 
Cluster results are now available for several such nuclei 
and we are interested in seeing how well the JDMRG 
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truncation strategy works on these nuclei and how well 
it compares with the Coupled Cluster method. 

Another issue of current interest concerns the possibil- 
ity of breaking up large single-particle orbitals, rather 
than adding them in a single stage. We have some 
thoughts on how this might be done without losing an- 
gular momentum conservation and we are now testing 
these ideas on the /7/2 orbital for the same nuclei we 
have already studied. If successful, we will then consider 
the application of these ideas to even larger shells, as will 
be critical for subsequent applications of the method to 
heavier nuclei, the ultimate goal of this project. 
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